A1 = readmatrix('LD_633_Pol.csv')

x = (A1(:,1)*2-17)*2*pi/360
y = abs(A1(:,2))
z = A1(:,3)
y= y./z

x1 = (A1(1:32,1)*2-17)*2*pi/360
y1 = abs(A1(1:32,2))
z1 = A1(1:32,3)
y1= y1./z1

x2 = (A1(32:78,1)*2-17)*2*pi/360
y2 = abs(A1(32:78,2))
z2 = A1(32:78,3)
y2= y2./z2

x3 = (A1(78:122,1)*2-17)*2*pi/360
y3 = abs(A1(78:122,2))
z3 = A1(78:122,3)
y3= y3./z3

x4 = (A1(122:167,1)*2-17)*2*pi/360
y4 = abs(A1(122:167,2))
z4 = A1(122:167,3)
y4= y4./z4

x5 = (A1(167:181,1)*2-17)*2*pi/360
y5 = abs(A1(167:181,2))
z5 = A1(167:181,3)
y5= y5./z5

ax = polaraxes
polarplot(x1,y1,LineWidth=2,Color=[0 0 0.8])
hold on
polarplot(x2,y2,LineWidth=2,Color=[0.8 0 0])
polarplot(x3,y3,LineWidth=2,Color=[0 0 0.8])
polarplot(x4,y4,LineWidth=2,Color=[0.8 0 0])
polarplot(x5,y5,LineWidth=2,Color=[0 0 0.8])
polarscatter(x,y,5,'MarkerEdgeColor',[0 0 0])
%rlim([0,1.1e-5])
ax.ThetaZeroLocation = 'top'
ax.ThetaDir = 'clockwise'
hold off

set(gcf,'Position',[0, 0, 200, 200])